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Abstract 

We consider operators of the form C = —L — V, where L is an elliptic operator 
and V is a singular potential, defined on a smooth bounded domain ft C lR n 
with Dirichlct boundary conditions. We allow the boundary of f2 to be made of 
various pieces of different codimension. We assume that C has a generalized first 
eigenfunction of which we know two sided estimates. Under these assumptions we 
prove optimal Sobolev inequalities for the operator £, we show that it generates 
an intrinsic ultracontractive semigroup and finally we derive a parabolic Harnack 
inequality up to the boundary as well as sharp heat kernel estimates. 

AMS Subject Classification: 35K65, 26D10 (35K20, 35B05.) 
Keywords: Singular heat equation, optimal Sobolev inequalities, generalized 
first eigenvalue, parabolic Harnack inequality, heat kernel estimates, degenerate 
elliptic operators. 



1 Introduction 

Let VL C M n , be a bounded domain and suppose V is an Lj oc (Q) potential for which 
we know the following L 2 estimate 

uGC °°(n) J n u 2 dx 
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One of the motivations of the present work is whether one can improve the above 
estimate to a Sobolev type estimate, involving, if possible, the critical Sobolev exponent. 
It is clear that to improve the previous estimate one needs more information concerning 
the potential V, besides (1.1). One additional piece of information that we are going 
to use, is the existence of a generalized eigenfunction <f>\ of problem (1.1) as well as 
sharp two sided estimates of <p\. Under this extra piece of information we are able to 
obtain sharp Sobolev type inequalities involving the critical Sobolev exponent. 

The knowledge of the asymptotic behaviour of <j)\ usually comes as a consequence 
of the maximum principle and the local character of V . We present such an argument, 
following the ideas of Brezis, Marcus and Shafrir [5] for the critical potential V{x) = 
\ in the Appendix. We should mention that in this case the asymptotics of <pi 

have already been derived by Davila and Dupaigne [12], [13]. The argument we present 
is much simpler and is based only on the maximum principle applied in the appropriate 
energy space. All the potentials we have in Section 4 as well as other potentials can be 
treated similarly. 

Presupposing the existence and asymptotic behaviour of the generalized eigenfunc- 
tion 4>i seems to be a natural assumption. It is E. B Davies [8] who put forward the 
idea of connecting the asymptotics of 4>i to the asymptotics of the Green function in 
the case of subcritical potentials V. In fact, he conjectured that knowing that the 
asymptotic behavior of 4>\ is like d{x) := dist(x, dVt) is actually equivalent to the exact 
two-sided Green function bounds. This conjecture was answered positively in [20] even 
for a larger class of potentials for which the generalized eigenfunction <pi behaves like 
d a (x) := dist Q (x, d£l), for some a > 1/2. 

The idea of obtaining heat kernel estimates of second order elliptic operators with 
singular potentials in terms of the generalized ground state is not new and besides [23] 
and [24] has been successfully exploited in [26] and [27]. 

Another motivation is the study of the corresponding parabolic problem, especially 
when the potential V is singular (see e.g. [1], [6] and [31]). In connection with this, we 
mention the work of Cabre and Martel [6] where the condition Ai > — oo is shown to 
be necessary and sufficient for the existence of a global positive weak solution. They 
also show that these solutions grow at most exponentially in time for any nonnegative 
initial data no G L 2 (VL). In fact, as far as the parabolic problem is concerned, condition 
(1.1) is practically equivalent to Ai > — oo due to a shift in the time variable. Under 
the extra assumption that the asymptotics of the generalized eigenfunction are known, 
we show that the corresponding Schrodinger operator generates a semigroup of integral 
operators for which we obtain parabolic Harnack inequality up to the boundary and 
precise heat kernel estimates. 

At this point we introduce some notations that we keep throughout the work. 
We assume that Q, C M n , n > 2, is a bounded domain, with a boundary <9S1 = 
U^ =1 Tfc, where = U^r^j is a finite union of m k distinct smooth C 2 boundaryless 
hypersurfaces T^j of codimension codimTfcj = k, where j = 1, . . . , k = 1, . . . , n— 1; 
in addition T n = {aq, . . . ,x mn }, T k j n = if k ^ /, or i ^ j and T k j n T n = 0, 
for k = 1, . . . , n — 1 and j = 1, . . . , m k . We also set d k {x) = dist(x, T^). For x £ Q 
we denote by d{x) the distance to the boundary dQ.. We clearly have that d(x) = 
min x£ Q{di(x), . . . , d n (x)}. Finally for the part of the boundary that is of codimension 
one we use the special notation diQ, = T\. 
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We are interested in the quadratic form 



L ( E 

il \i,j=l 



Q[y] = / J E a ^^§^§^ ~ Vu * \ dx > ue ^(n), (1-2) 



where V G L\ oc {VL) and ajj(x), i, j = 1, . . . , n is a measurable symmetric uniformly 
elliptic matrix, that is 

n 

Co|£| 2 < E < Co~W> e e F (1.3) 

for some Cq > 0. We also assume the L 2 estimate 

0<A i: = inf (1.4) 

and that to Ai, there corresponds a generalized eigenfunction </>i of (1.4). More pre- 
cisely, we assume that 4>\ G Hj oc (VL) n L^ c (r2) and that 



X e a ^4t-^ dx = L {v + x ^ Mdx > ^^(A). (i.5) 

In addition we assume that we have two sided estimates on <pi of the form <fi± ~ 
d? 1 ...d£", that is 

cid? 1 ^)...^"^) <^i(x) ^cad? 1 ^)...^"^), (1.6) 

for any x G 0, for two positive constants ci,C2 and for suitable exponents a\, . . . ,a n . 
Appropriate conditions on the exponents ai will be formulated below. 

Our first result concerns the following improved Sobolev type inequality. 

Theorem 1.1 (Optimal Sobolev type inequality) For V € Lf oc (Q) we assume 
that (1-4) holds and in addition there exists a ground state <fii G Hi oc (Q) n Lf£ c (Q) 
satisfying (1.6) for 

afc> -^- (n - fc) 2^rV fe = 1 '-' n > C 1 - 7 ) 
where 2 < q < if n > 3 and g > 2 if n = 2. We then have that 

0<C(n,a 1 ,...,a n ,q)= inf 9M (1. 8 ) 

«6C -(Q) / „(n-2)-2n \J 

In particular when n > 3 and g = we ^cwe i^ai 

< C(n, ai ,...,a n ) = inf Q ^ 



(ln\ u \ n - 2dx ) 
if 

k — 2 (n-k) 

dk> 7, 777 fv. fc = l,...,n. 

2 2(n — 1) 
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We note that the condition (1.7) is optimal. In the extreme cases q > 2 and 
a n = or q = 2 and at = — ^-^ we have different improved inequalities, see 

Theorems 2.4 and 2.5 respectively. 

Using the above Sobolev type inequality we will now proceed to the study of the 
corresponding parabolic problem, that is 

f % = -£u:=EZ j =i 1 §r i (a ij (x)J^)+V(x) U in (O.ooJxO, 

I u(x,t) = on (0,oo)xd 1 n, (1.9) 

[ u(x,0) = 1*0(2;) on . 

Our first result is the following Harnack inequality. 

Theorem 1.2 (Parabolic Harnack inequality up to the boundary) For V G 

Lj 1 c (f2) we assume (1-4) and that (1.6) holds for some > for k = 1, ■ ■ ■ , re. 

Then for C as in (1.9) there exist positive constants Cu and R = R(Q) such that for 
x £ £1, < r < R and for any positive solution u(y, t) of 

— = -Cu in {B{x, r)nn} x (0, r 2 ) , 
the following estimate holds true 

n 

sup (2/t)e{B(:Ci§)nn}x( ^ ^ ) u(y,t)l[d- ai (y) < 

i=l 

n 

< C u ess mf (?/it)g { B(;Cii)nn } x( 3 r2ir2) u(y,t) J[ d~ ai (y). 

i=i 

Here B(x, r) denotes roughly speaking an n dimensional cube centered at x and 
having size r, see Definition 3.1 for details. 

Theorem 1.2 states a parabolic Harnack inequality up to the boundary for the ratio 
of any positive local solution to the Cauchy-Dirichlet problem and the generalized 
eigenfunction (f>\. We note that a\ > 1/2, therefore <p\ is zero on the boundary d\VL. In 
particular it implies that any two nonnegative solutions vanishing on <9iO must vanish 
at the same rate. It is clear then that such a normalization is necessary. In fact the 
natural quantity is v = u/(f>\ and it is for this function that we prove the Harnack 
inequality. For the definition of solution for the function v we refer to Definition 3.8, 
where however the appropriate weight is <f>\. 

Alternatively, one could define local weak solutions of (1.9) using a suitable local 
energetic space obtained via the quadratic form (1.2). For an example of globally 
defined energetic solutions see [31]. 

Our result in the case a.\ = 1, at = for k = 2, ■■■,n is basically the local 
comparison principle of [14]. We also note that the restriction on in Theorem 1.2 is 
sharp. 

In what follows we denote by h the integral kernel of the 1? semigroup associated 
to the elliptic operator C as defined in (1.9), that is 

u(x,t) := / h(t,x,y)u (y)dy . 
Jq 
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The existence of h(t, x, y) is proved in Proposition 2.8 and it is a consequence of our 
Theorem 1.1. 

As usual, from the parabolic Harnack inequality one can obtain sharp heat kernel 
estimates, as explained by Grigoryan and Saloff Coste, see [23], [24], [30]. In particular 
we have 

Theorem 1.3 (Sharp heat kernel estimates) For V G Lj oc (Q) we assume (1-4) 
and that (1.6) holds for some at > — f or k = 1, ■ ■ ■ ,n. Then there exist positive 
constants C\,Ci, with C\ < Ci, and T > depending on such that 

for all x,y £ fl and < t < T, whereas 

n n 

d J] d? (x) df (y) e- Alt < h(t , x, y) < C 2 J] df (x) df (y) 

i=i i=i 

for all x,y £ Q and t >T. 

Due to a shift in the time variable, Theorems 1.2 and 1.3 remain valid if we replace 
assumption (1.4) with the condition Ai > — oo. For the corresponding statement of 
Theorem 1.1 under the condition Ai > — oo we refer to Theorem 2.1. 

Although we present here only heat kernel estimates one can integrate in time to 
obtain the corresponding Green function estimates provided that Ai > 0. 

It is clear that the asymptotics of <f)\ affect both the parabolic Harnack inequality 
and the heat kernel estimates, see Theorems 1.2 and 1.3. At the same time it seems 
that the Sobolev inequality is independent of the a^s, in the sense that the exponents 
afc do not appear in the ratio (1.8). We note however that there are critical cases where 
relation (1.8) fails and different Sobolev inequalities hold true. For instance if a n = 
— ^—^ estimate (1.8) is no longer true; instead, the optimal Sobolev inequality involves 
a logarithmic correction see inequality (2.23) in Theorem 2.4. For other examples see 
Theorem A' in [22]. 

We note that instead of the uniform ellipticity condition (1.3) which we assume 
throughout this work, our method can also treat degenerate operators for which the 
following condition holds 

Cow(x)\Z\ 2 < J2 aij(xMj < C^wixM 2 , £ G M n , (1.10) 

where w(x) is like a power of the distance function. 

Finally we should mention that this work complements and extends our previous 
work [20]. There we studied the cases where the potential V(x) is either for a 

general bounded domain Q, or else , .. ■>) — ^=rr for a convex bounded domain Q. In the 
second case the convexity was used in an essential way. Here, even in these two cases, 
we improve our results in the first case by allowing potentials involving distances to 
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more than one point and in the second case by removing the convexity assumption (see 
Section 4). 

The article is organized as follows. In Section 2 we establish a Sobolev type in- 
equality, thus giving the proof of Theorem 1.1, starting from the 1? estimate and using 
the behavior of the generalized eigenfunction (pi . In Section 3 we study the associated 
Cauchy-Dirichlet problem and prove a parabolic Harnack inequality up to the boundary 
as well as sharp two sided estimates on the corresponding heat kernel. In particular we 
provide the proof of Theorems 1.2 and 1.3. Finally, in Section 4 we give some examples 
of concrete Schrodinger operators with singular potentials for which an L 2 inequality 
holds true and the behavior of (pi is known. In all these examples the results of the 
present work apply. 



2 From the L 2 estimate to Sobolev type inequalities 

In this section starting from the 1? estimate (1.4) we will prove various Sobolev inequal- 
ities involving optimal exponents. In particular we will prove a weighted logarithmic 
Sobolev inequality that will be crucial in establishing the intrinsic ultracontractivity of 
the semigroup associated with the operator C defined in (1.9). 
For 5 small enough we set 

rf = { x g n, s.t. dist(x, r fc ) < 5} and (r£) c = n \ r{ . 

As a consequence of the assumptions on the domain Q we made in the introduction, 
we have that for 5 small enough r<^. n rf- = if A; ^ I, or i ^ j and T s k j n T 5 n = , 

for k = 1, . . . ,n - 1 and j = 1, . . . ,m k . We note that n = (ug =i rjQ U (u£ =i r£) C = 

(ug = iit)u(nj =1 (it) c ). 

We are now ready to give the proof of Theorem 1.1. Indeed we will more generally 
assume instead of (1.4) the following I? estimate 



-oo<Ai:= inf (2.1) 



Then we prove: 



Theorem 2.1 (Optimal Sobolev type inequality) For V G Lf oc (Q) we assume 
that (2.1) holds and in addition there exists a ground state (pi G Hl oc (Q) n Lf£ c (Q) 
satisfying (1.6) for 

ak> -±^- in - k )j£L., k = l,...,n (2.2) 

where 2 < q < if n ^ 3 and q > 2 if n = 2. We then have that for every A > 

,n „ Q\u] + (A — Ai) f u 2 dx , n „. 

< C{n,ai,...,a n ,q,X) = inf ^ -. (2.3) 

ueC™(Q) / „(„-2)-2n \q 

[Jnd 2 \u\idx\ 
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In particular when n > 3 and q = we have that for every A > 0, 

Q[u] + (A - Ai) f n u 2 dx 



if 



< C(f2, ai, . . . ,a n , X) = inf 



k — 2 (n-k) 

&k > X 7w tt, fc=l,...,n. 

2 2(n — 1) 



In the case Ai > one can take A = Ai, thus proving Theorem 1.1. 
Proof of Theorem 2.1: It is a consequence of the following estimate for any v € Cq°(Q) 



J n ^yp a n v *i v xi + At,2 j dx - c (j n ^ q{n * 2n \ v \ 9dx ) q . ( 2 - 4 ) 

with C = C(qi, . . . , a n , f2, A) > 0. For convenience we write v Xi instead of Let 
us accept (2.4) and give the proof of the Theorem. Clearly (2.4) is valid not only for 
smooth functions but also for functions in the completion of C^°(Q) under the norm 
defined by 

\\v\\ Hli := ^Jj 2 (v 2 + \Vv\ 2 )dx^ . (2.5) 
In particular we can take v = with u € Cq°(Q), in which case we get 

A f , „f / (<t>l)x s A , ^ f (MxMxj 2 , 

> / aaUx-Ux-dx — I > / aijuu x . — ■ — i ax+ > / a;,- -u dx 

i,3 h3 h3 1 

2 

u 2 dx>C(J d s | u |9dxJ . (2.6) 

2 

On the other hand, by standard approximation arguments, (1.5) is valid also for tp = 
with u £ Cq°($1). For this choice of the test function, we get that 



2 > / aijuu x - -dx — } / ajj 5 n dx = (V + Ai)u dx. (2.7) 



<j)l x ^ 

1,3 h3 

Combining (2.6) and (2.7) we conclude that for any u 6 Co°(f2), 



/■ „ / f q(n-2)-2n \ q 

Q[u] + (A - Ai) / u 2 dx >C[ d 5 , 

which is the same as (2.3). 

It remains to prove (2.4). Because of the ellipticity condition (1.3), estimate (2.4) 
follows from 

4>l (\Vv\ 2 + v 2 ) ris > C ^Id 1 ^ 2 ^ \v\ q dx^j " ,veC%°(n). (2.8) 

In view of (1.6) we may replace 0i in (2.8) by d" 1 . . . d^ n ■ Estimate (2.8) then is true 
in n£ =1 (r£) c . This is a consequence of the standard Sobolev imbedding of functions in 
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# 1 ( n fc=i( r k) c ) in Lq ( n k=i( T k) c ) and the fact that 6 < d k( x ) < D k ■= sup xen d k (x) < 
oo, k = 1, . . . , n. 

We therefore need to prove that estimate (2.8) is true when we replace by UjJ =i r|. 
As a matter of fact it is enough to prove that for any k = 1, . . . , n, 



J d 2 k ak (\Vv\ 2 + v 2 )dx > C \J d% kq \v\ q dxj \ v £ C °°(ft), 

where f3 k = a k — 1 + ^^n. The validity of this estimate is given in the next main 
Lemma. 

Lemma 2.2 Let $7 C M n , n > 2, be a bounded domain. Suppose that T k C d£l is a 
smooth boundaryless hypersurface of codimension k, k = 1, . . . , n — 1. When k = n we 
take T n to be a point. We also assume that 

f3 k = a k -l + (2.9) 

where 2 < q < if n > 3, q > 2 if n = 2. Then, for any k = 1, . . . , n, there exists a 
C = C(a k , 5,k) > such that for all 5 > and all v £ Co°(J7) there holds 

J d 2 k ak (\Vv\ 2 + v 2 )dx > C\\d^v\\ 2 Lq{Tiy (2.10) 

fc 

provided that 

k — 2 a — 2 

dk + {n-k)^- r. (2.11) 

^ 2 v 2(<? + 2) v ; 

Proof: Let us fix a k = 1, . . . , n. We will initially establish the result for 5 small. For 
simplicity we write d instead of d k . From Lemma 4.2 [19] (see also [18]), we have that 
if 

n — 1 

KQ< , b = a-l + ^——n, and a^l-k, (2.12) 

n — 1 Q 

then, for S small there exists a C > such that there holds 

C\\d h w\\ LQ(V 5,< [ d a \Vw\dx+ [ d a \w\dS x , weC^°{n\T k ). (2.13) 

We apply (2.13) to the function w = \v\ s , s = Also, for a k , (3 k , and q as in (2.9), 
we set 

Q = qs , b = [3 k s, a = 6 + 1 — n = — + a fe . 

It is easy to check that a, 6, Q thus defined satisfy (2.12). As far as the condition 
a 1 — k is concerned, when written in terms of a k , q, k and n, it is equivalent to 

a k T 7, \ n ~ k ): 



2 y 2^ + 2)' 
which is precisely (2.11). From (2.13) we have 



nv >:!; : , = c\\d b \v\ s \\ LQ(r s k) < s f d a \v\ s ~ 1 \v v \dx+ f 6 d a \v\ s ds x , (2.14) 
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for some positive constant C. Using Holder's inequality in the gradient term of the 
right hand side we get 



/ d>| a-1 |Vu|da; = / d ak \Vv\ d—\v\*dx 
M Jrl 



< \\d ak \Vv\\\ L2{Ti) \\d^v\\l q(Tl) 

< ^KMVwIII^+ell^wll^. (2.15) 
Hence, from (2.14) and (2.15) we arrive at 

(C - es)\\d^v0 Ti) < sc £ \\d a * {VvlW 1 ^ + J d a \v\ s dS x . (2.16) 

k 

To continue we will estimate the trace term in (2.16). Using Holder's inequality we 
have that 



/ d a \v 
JdTl 



(n-2)(2 + g) 

2+9 / f 2(n-l) \ 4(n-l) 

\ s dS r = 5 a \v\ — dx<c5 a \v\^=^dx) 



9T k \ Jdr l 

Af>k-<*k)<! 1+1 

L-s=3-(ar«) 



By the trace imbedding [3], we have that for u G If 1 (r|' 



(2.17) 



2 ^ /-</„ L\iir7».ll2 , )irii..ii2 



W 1 1 2(n-l) 



C(n,A;)||V«||i 2(ri) + M||«||^ (ri) , (2.18) 



l n-a (ar*) k 
where M = M(n, r£). Applying this to u = <i afc i; we get 

||cT fc v|| 2 2(n-i) < C 2 / d 2ak v 2 dx + C 2 I d 2ak \Vv\ 2 dx, (2.19) 

with C 2 = C 2 (afc, (5, fc, n). 

Indeed after some elementary calculations from (2.18) we get for any 9 > 1 

IK fc+e H| 2 2 ( n-i) < 2C(n,k)(a k + e) 2 f d 2 ^ k+e ^' 2 v 2 dx 

L^^(drl) Jrl 

+2C(n,k) [ d 2 ^ ak+ ^\Vv\ 2 dx + M f d 2( - ak+e ^v 2 dx 

< {2C{n,k)(a k + efS 29 ' 2 + M5 2e ) [ d 2ak v 2 dx 

H 

+2C(n,k)5 2e ( d 2ak \Vv\ 2 dx. 

Whence, 

N afc ^|| 2 2 ( »-i) < (2C(n,k)(a k +9) 2 5- 2 +M) f d 2ak v 2 dx+2C(n,k) f d 2ak \Vv\ 2 dx, 
L^=^{dr{) Jrl Jr s k 

(2.20) 
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and therefore (2.19). We note that (2.19) is valid even for nonpositive values of a k . 

Combining (2.16), (2.17), (2.19) and then raising to the power we easily con- 
clude (2.10) for 5 small. 

The general case follows by noticing that outside for small 5 the corresponding 
estimate comes from the standard Sobolev embeddings and the fact that 5 < d k {x) < 
D k := sup xen d k (x) < oo, k = 1, . . . , n. 

This completes the proof of the Lemma as well as of Theorem 2.1. 

□ 

We note that when k = n condition (2.11) reads a n ^ —^^ L . It turns out that the 
analogue of the estimate (2.10) in case k = n and a n = — involves logarithmic 
corrections. More precisely we have: 

Lemma 2.3 Let £1 C M n , n > 3, be a bounded domain and T n = {xq} C d£l be a 
point, such that for some 5q small Tf° = Bs () (xo) \ {xq} C £1. We also assume that 

n-2 n 2n n n-2 i q-2 

«n = x—; 2<q< -, n = 1 + — — n. 

2 n — 2 2 2q 

Then, there exists a C = C(Q,5, n) > such that for all 5 > and all v £ Cg°(f2) 
there holds 



d 2 - n (\Vv\ 2 + v 2 )dx > C\\d^X-^v\\ 2 Lq{vi) , (2.21) 



where X = X{^£L), with X{t) = (1-lnt) -1 , < t < 1, and D n := sup xen d n (x) < oo. 

Proof: As in the previous Lemma it is enough to give the proof for 5 small. We may 
assume that xo = 0, hence, d n {x) = \x\. Also, for simplicity we suppose that 5 = 1. 
Then we recall the following result for any w G Cg^f^) the following estimate holds 

2 

\x\ 2 ~ n \Vw\ 2 dx>c(^J B \xf nq X 1+ 2\w\ q dxY =c\\d l3n X^w\\ 2 Lq{B2) . (2.22) 

This is Lemma 3.2 [22] in the case q = an d Proposition 6.2 [2] in the case where 
2 < q < Given a function v G Cq°(Q), we extent it from B\ to the function v 

supported in P>2 so that 

\\v\\b 2 < C\\v\\ Bl , 

holds for a positive constant depending only on n, where we denote by \\v\\% := 
f Bi \x\ 2 ~ n (\ Vf | 2 + \v \ 2 )dx (note that away from the origin v is an H 1 function). We 
next apply (2.22) to v and the result follows easily. We note that one cannot take a 
smaller exponent of the logarithmic term X. 

a 

As a consequence of the above Lemma we have: 

Theorem 2.4 For V £ L} oc (tt) we assume that (2.1) holds and in addition there exists 
a ground state <f>\ G Hl oc (Q) n Lf£ c (p.) satisfying (1.6) for k = 1, . . . , n — 1; n > 3 

k-2 . ,, q-2 n-2 n 2n 
a k> 7, (rx — A;) — — — , a n = — , 2<q< 



2^ + 2)' n 2 ' H ~n 
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We then have that for every A > 

0< <?<*«, «.,,,*)- inf OH-KA-AOJn^ ^ (MJ) 

where X = X(^ff-), with X(t) = (1 - Int) -1 , < t < 1 and L>„ := sup^n d„(x) < oo. 
In particular by choosing q = we /iaue i/iai /or every A > 

o ^,r^ ^\ • r <9[«] + (A - Al) / n ti 2 dx 

0<C(n,ai,...,a n ,A)= mf -^U ^ Jm n _ 

U€C%°(n) / 2(n-l) 2n X ~ 

J n X -2 \u\»-*dx 



Proof: The proof is quite similar to the proof of Theorem 2.1, where in the place of 
Lemma 2.2 one uses Lemma 2.3 for k = n. We omit further details. 

□ 

Concerning the limit case q = 2 and = — we have the following 

Theorem 2.5 For V G -£j oc (ft) w/e assume that (2.1) holds and in addition there exists 
a ground state <f>i G Hi oc (£l) n Lf£ c (£l) satisfying (1.6) for > —^r, k = 1, . . . ,n; 
n > 2. We then have that for every A > 

n/^o u ■ t QH + {\-\i)!n u2dx 
< C(S2, «i, . . . , a n , A) = mf 5 — , 

w/iere X = X{^^-), with X(t) = (1 — lni)^ 1 , < t < 1, and D := sup xe Qd(x) < oo. 

Proof: The proof is quite similar to the proof of Theorem 2.1, where in the place of 
Lemma 2.2 one uses Lemma 2.6 below. We omit further details. 

□ 

Lemma 2.6 Let Q C M n , n > 2, be a bounded domain. Suppose that T/, C d£l is a 
smooth boundaryless hypersurface of codimension k, k = 1, . . . ,n — 1. When k = n we 
take r n to be a point. Then, for any k = 1, . . . ,n, there exists a C > such that for 
all 5 > and all v G Cq°(Q) there holds 

I d k {k ~ 2 \\Vv\ 2 + v 2 )dx >C [ X 2 d k ~ k v 2 dx, 

where X = X( dk ^ ), with X{t) = (1 — hit) -1 , < t < 1, and Dk ■= sup xen dk(x) < oo. 
Proof: This is proved using similar ideas as in Lemma 2.3. We omit further details. 

□ 

As a consequence of Theorems 2.1 and 2.4 we obtain 

Theorem 2.7 (Weighted log Sobolev) For V G L} oc (tt) we assume that (2.1) holds 
and in addition there exists a ground state 4>\ G Hl oc (£l) n Lf£ c (Q) satisfying (1.6) for 
n > 2 with 

k — 2 n-k . 1 ^ n-2 , nni \ 

"2(^1)' fc = l,...,n-l, a«>-— (2.24) 
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Let 

A := max{ai, (X2, ■ ■ ■ , a n , 0}. 

Then, there exists a positive constant K such that for any e positive and for any u G 
Cg°(fi) there holds 

i uha ( mi J""'... as- ) dx - eQ[u] + { K ~ n - ± r L ln£ ) Ml < 2 - 25 > 

/iere H^H 2 , = / n |ii| 2 (i:r. 

Proof: At first we will show that in each r|, the following estimate holds: 

J^\ v ^^M^jdx<e <t>\(\Vv\ 2 + v 2 )dx +^K- U + 4 2 "^ lne^j ||v||£ fc , (2.26) 

where := maxja^., 0} and |Mlr fc := Jp* <f>\V 2 dx. To this end, let us assume first that 
u; is normalized so that for d/j, = <\>\w 2 dx one has / r * dfi = ||to|| r = 1. Then, for q > 2, 
using Jensen's inequality, we have 



2 



I d>W MM dx = ^~ 2 I ln\ W r 2 d, < ^^In I <plH q dx\ " . (2.27) 
To continue, we will use the estimate 

2 

( ^ tfM'd* 9 <C J s <p\{\Vw\ 2 + u> 2 )cfe. (2.28) 

In case fc = 1, ... ,n — 1 or k = n and a ra > — (2.28) is a direct consequence of 
(2.10), provided that d 2 k ak < cdf k , that is, 2a k > qf3 k . In view of the definition of f3 k 
(see (2.9)), the requirement 2a k > q(3 k is equivalent to 

q(n-2 + 2a k ) < 2(n + 2a k ). 

We note that when n > 3 if a k < 0, then we can choose q = whereas if a k > 0, 

then the maximum q one can choose is q = ^j^2+2a k ) ^ n^2' ^he same choice of q is 
feasible when n = 2 and a±, a>2 > 0. Hence, in any case one takes 



2(n + 2a+_^_ 



(n - 2 + 2a; 



(2.29) 



A; 



On the other hand, in case k = n and a n = — ^-^y 2 -, estimate (2.28) is a direct conse- 
quence of (2.21) with q = if n > 3; indeed, in this case one has q(3 n = —n and 

clearly <i 2_ra < cd~ n xi +1 . In particular, in all cases the choice of q is given by (2.29). 
^From (2.27) and (2.28) we get that 

/ 0jw 2 ln\w\dx < , q A n [C [ <j>\{\ Vw\ 2 + w 2 )dx j . 

Jv{ 2( 9 -2) y / 
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Using the fact that ln# < sd — lne for all 6, e positive we get that there exists a K > 
such that for any e > 

/ <t){w 2 \n\w\dx <e ( 4>l(\Vw\ 2 + w 2 )dx + K - - g . lne. (2.30) 
Jv{ Jv{ 2(g-2) 

Because of (2.29) we have that 

q _ n + 2a~l 
2(9 - 2) " 4 • 

On the other hand given any f G Cg°(f2) we apply (2.30) to w = to conclude 

(2.26). 

We next consider aw £ Cq°(Q) normalized by \\w\Iq := J n (p 2 w 2 dx = 1. Applying 
(2.26) to this w; we get 

/ </> 2 u; 2 ln|tt>|<ix - / 4>\w 2 h\(\\w\\r k )dx < e / 2 (| Vu>| 2 + w 2 )c?x 
7r 1 M 

n+ 2a+ V 2 
+ ^ ^-IneJHI^. 

Since ||u> ||r fe < 1 and therefore ln(||iy||r fc ) < 0, we have in particular that 
/ (t)\w 2 \n\w\dx < e [ 2 (|Vu>| 2 + w 2 )dx + ( K 

H J n V 

Summing over all Tf. we get 

/ (ftwhn^dx < e f (j) 2 (\Vw\ 2 + w 2 )dx +(k~ - + ? A lne) IMl?,^. (2.31) 
Jur* Jurt \ 4 / 



F^ lne Hlr fc - 



On the other hand on f2 \ ur| we have that <pi ~ C and using the standard log Sobolev 
inequality we easily arrive at 

/ 4>\w 2 \n\w\dx < e [ 2 (| Vw\ 2 + w 2 )dx + ( K - - t lne J |M| 2 ur ^, (2.32) 
J(ur^) c i(ur*) c V 4 / v fc; 

when n = 2 (2.32) holds true for any v > 2 in place of n. Combining (2.31) and (2.32) 
we get that for \\w\\q = 1, there holds 

J^4>\w 2 Mw\dx <e J^\{\Vw\ 2 + w 2 )dx+ (k- n+ ^ A in^ . ( 2 .33) 

For a general u G Co°(r2) we apply (2.33) to w = to get 

/ ^yin ( H ] dx < £ f 2 (| Vw|2 + u2)da . + ^ _ \ 2 _ (2 34) 

./n \IM|n/ ./n V 4 / 

Taking u = ^u, (2.25) follows. 

□ 

The above logarithmic Sobolev inequality is the main ingredient in establishing the 
intrinsic ultracontractivity of the semigroup generated by the operator C defined in 
(1.9). More precisely we have 
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Proposition 2.8 Let V G L\ oc (Vl). We assume that (2.1) and (1.6) hold for some 
a k > —^-^ — 2(n-i) > ^ = • • • ' n ~ ^> a ™ — ~~ *T^ - Then the operator C defined in 
(1.9) gives rise to an intrinsic ultracontr active semigroup in L 2 (Q), whose heat kernel 
h(t, x, y) satisfies 

h(t, x, y) < c £Ml^!(g^hM e -A lt for anyt>0,x,ye{l; (2.35) 

t 2 

where A := max{a\, ct2, ■ ■ ■ , a n , 0}. 

Proof: This is quite similar to Theorem 3.4 in [20] for this reason we only sketch it. 
We change variables by 

v{x,t) :=u{x,t)/fa(x), (2.36) 
then if u solves problem (1.9) the function v satisfies 

f ft = -^i«:=^E? J -=ilfe(^ a y(^)-Ai« m (0,oo)xn, 

{ v(x,t)=0 on (0,oo)xd 1 n, ( 2 - 37 ) 

[ v(x, 0) = wo(^) on f2 , 

with ■uo(x) := tio(^)^^ 1 (x). 

We note that the elliptic operator — X\ is defined in the domain Dl^C^ — Ai) := 
{v e Hq(£1;4>i) : Cfa — Ai € L 2 (Q, <j>\(y)dy)} , where Hq(Q,;<PI) denotes the closure 
of Cq°(Q) functions with respect to the norm (2.5). To this elliptic operator it is 
naturally associated a bilinear symmetric form which is a Dirichlet form. Then Lemma 
1.3.4 together with Theorems 1.3.2 and 1.3.3 in [9] imply that the elliptic operator 
£>4>\ ~ -^l generates an analytic semigroup, e _ ^i _Al ^*, which is positivity preserving 
and contractive in L P (Q, 4>\dx) for any 1 < p < oo. 

^From the weighted logarithmic Sobolev inequality (2.34), we deduce the corre- 
sponding L p logarithmic Sobolev inequality for any p > 2; to this end it is enough to 

P 

apply (2.34) to the function v := for any smooth w. Using Theorem 2.2.7 in [9] - 
as it is used in Corollary 2.2.8 of [9] - we obtain that e - ^! - * 1 )* is an ultracontrac- 
tive semigroup. As a consequence the semigroup e~ c<t >i t has a heat kernel hfa, which 
satisfies the following uniform upper bound 

C 

V (t, x, y) < -^xe~ Xlt for any t > 0, x, y G n . (2.38) 

t 2 

Clearly the heat kernel upper bound (2.35) associated to the operator C follows from 
(2.38), (1.6) and the fact that 

h(t, x, y) = 0i(x)0i(y)/t^ (t, x, y) , (2.39) 

which is an immediate consequence of the change of variables (2.36). We omit further 
details. 

□ 
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3 Harnack inequalities and sharp heat kernel estimates 



In this section we prove a parabolic Harnack inequality up to the boundary for the 
operator defined in (2.37), and we deduce from it the corresponding heat kernel 
estimates as well as the proofs of Theorems 1.2 and 1.3 in the Introduction. We use 
Moser iteration technique, as adapted in [20] for bounded domains fl. To this end we 
will prove four basic estimates. Namely, a sharp volume estimate, a local weighted 
Poincare inequality, a local weighted Moser inequality and a density theorem. 

We will use the following local representation of any smooth boundaryless hypersur- 
face Tfcj of codimension k = 1, . . . , n — 1, for any fixed j = 1, ■ ■ ■ , mt, which is of course 
Lipschitz. That is, we suppose there exists a finite number ./V (depending on both k 
and j) of coordinate systems (yi,Zi), yi = (yi,i, ■ ' ' ,Vi,(n-k)) and *i = (zi,i, ■ ■ ■ , *i,k), 
for i = 1,---,7V, and the same number of functions a« = dj(yi) : M n ~ k — > M k , 
(aj = (aj,...,a k )) defined on the closures of the (n — k) dimensional cubes Aj := 
{Vi '■ \Ui,i\ < (3 for I = 1, ■ ■ ■ , n — k}, i G {1, • • • , N} so that for each point x G T^j 
there is at least one i such that x = (yi,CLi(yi))- The functions a« satisfy the Lipschitz 
condition on Aj with a constant L > that is 

\ai{yi) - ai(yi)\ M k < L\y { - yi\ M n-k 

for any yi,yi G Aj. We note \y\mk is the Euclidean norm in M k . Moreover, there exists 
a positive number (3 < 1, called the localization constant of F^j and S7, such that the 
set Bi defined for any % G {1, ■ • • , A^} by the relation 

Bi = {(yi, Zi) : yi G \, a\( yi ) - (3 < z iyl < a 1 ^) + 0} , 

satisfies 

Ui = Bi n n = {( yi , : Vi G A { , a\(yi) - (3 < z^ < a\(yi)} if k = 1 , 

or 

f/i = B, t nn = Bi = {(yi,Zi) : Vi G A, t , a\( yi )-{3 < z^ < a\( yi )+(3} if k = 2,...,n-l , 

and Ti = B{ (1 dfl = {(yi, z^ : yi e Ai, Zi = (ii(yi)}. Finally, let us observe that for any 
y G Ui one has 

(1 + L^la^yi) - Zi\ M k < d k (yi,Zi) < \ai(yi) - Zi\ Mk ; 

see Corollary 4.8 in [25]. 

We fix a constant 7 G (1,2) and we define the "balls" we will use in Moser iteration 
technique. Roughly speaking they will be Euclidean balls if they stay away from the 
boundary and they will be n dimensional "deformed cubes" , following the geometry of 
the boundary, if they are close enough to the boundary or if they intersect it. More 
precisely we have 

Definition 3.1 (i) For any x G CI and for any < r < (3 we define the "ball" centered 
at x and having radius r as follows. B(x,r) = B(x,r) the Euclidean ball centered at x 
and having radius r if d(x) > jr (thus dk(x) > ^r for any k = 1, ■ ■ • , n) or if k = n, 
while 

B(x, r) = {(y { , z^ : \y { - x'\ R n-k < r, 
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a\(yi) - r - d k (x) < z it i < a\(yi) + r - d k {x) for any I = 1, ■ ■ ■ , k} 

if k = 1, . . . , n — 1 and d k {x) < jr where i G {1, • • • , N} is uniquely defined by the 
point x £ Tfc such that \x — x|jRn = d k {x), that is by the projection of the center x onto 
Tk C dCl, and x' denotes the first n — k coordinates of the point x in the i-orthonormal 
coordinate system, (ii) We also define the volume of the "ball" centered at x and having 
radius r by 

,. n 

V(x,r):= / \{dl a \y)dy. 



We first derive a sharp volume estimate, which plays a fundamental role in getting 
the sharp dependence of the heat kernel on x, y and t. 

Lemma 3.2 (Sharp volume estimate) Let n > 2 and a& > — | for k = 1, . . . ,n. 

Then, there exist positive constants c\,C2 and ro such that for any i£!! and < r < r^, 
we have 

n n 

ci I] + r ) 2akrH ^ V ( x > O < c 2 I] ( d k(x) + r) 2ak r n . (3.1) 

k=l k=l 

Proof of Lemma 3.2: Let us first consider the case where d(x) > jr, whence dk{x) > jr 
for any k = 1, ...,n. Then B(x,r) = B(x,r) C tt. Due to the fact that for any 
y £ B(x, r) and any k = 1, ■ ■ ■ , n we have 

(■^p) < 4(x) - r < 4(2/) < dfc(x) + r < (^^) 40*0 (3.2) 

we easily get 

n n 

V(x,r)~ r n J] 4 Qfc (x) ~ r n [] (4 (x) + r) 2a * , 
fc=i fe=i 

this proves the claim. 

Let us now consider the case where d(x) < jr. We claim that in this case there 
exists exactly one k = 1, ■ ■ ■ ,n such that dk{x) < jr. This is due to the assumption 
that for some 5 small enough r| ■ PI = for any k ^ I and i ^ j and r| ■ PI = 

for any = 1, • ■ ■ , n — 1 and j = 1, ■ • • , m^, since we may suppose that r < | (take 
ro := min{/3, |}, (5 being the localization constant of T^j and Q). Whence if dk(x) < jr 
then x £T s k and <i, (x) > 5 > jr for any j / fe, thus from (3.2) for any y G S(x, r) we 
have (y) ~ dj (x) ; as a consequence 

nc,r)~ I] (dj(x)+r) 2 ^ f df k {y)dy , 

j=ij& JB(x,r)nn 

Hence the claim will follow as soon as we prove that 

/ dl ak (y)dy~r n+2a *. (3.3) 

JB(x,r)nn 

Arguing as in (3.2) we have that dk{y) < (1 + j)r for any y G B(x,r). Moreover one 
has that if k ^ n, d}~(y) > r(j — 1) on a set of measure r n . Indeed 

/ dy= dz u dyi = 

JB(x,r)nnn{d k (y)>ri~f-l)} J\yi-x'\ Rk <r Ja\( yi )-r-d k (x) 
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= (2r - 7 r + d k {x)) k r n - k > (2 - 7 )V, 

and (3.3) follows. In the limit case k = n see [29] for any a n £ (— ^,0] (see also [28] 
and Lemma 2.3 in [20]). We note in fact that the same proof works for any a n > — ^. 

□ 

^From this one can easily deduce the doubling property: 

Corollary 3.3 (Doubling property) Let n > 2 and a k > — § for k = l,...,n. 

Then, there exist positive constants Cp> and r$ such that for any x £ Q and < r < ro, 
we have 

V(x,2r) < C D V(x,r) . 



Our next result reads: 

Theorem 3.4 (Local weighted Poincare inequality) Let n > 2 and a\ > 0, 

a k > — § for k = 2, ■ ■ ■ , n. Then, there exist positive constants Cp and ro such that for 
any iGS! and < r < ro, we have for all f £ C 1 (B(x, r) n fi) 

„ n „ n 

™l U d l ak (y)\f(y)-^dy<c P r 2 U d T k (y)\^f\ 2 dy ■ (3-4) 

£eK JB(x,r)nn ^ JB(x,r)nn ^ 

We note that our weight is not necessarily in the Muckenhoupt class A 2 . 

Proof: Let us first consider the case where d(x) > 7 r. Then B(x,r) = B(x,r) C SI 
and for any y £ B(x,r) and any k = 1, ■ ■ ■ ,ra we have d k (y) ~ dk(x), as in estimate 
(3.2). Thus in this case (3.4) follows from the standard Poincare inequality: 

inf/ \f{y)-i\ 2 dy<C P r 2 f \Vf\ 2 dy, f€C\Bfar)). 

e,^JRJB(x,r) JB(x,r) 

Let us now consider the case where d k {x) < jr, for some k = 1, ■ ■ ■ , n. Then arguing 
as in Lemma 3.2 it is enough to prove the following for any / £ C 1 (B(x,r) n SI) and 
any k = 1, ■ ■ ■ , n 

inf / \f(y) - t\ 2 d 2 k ak {y)dy < C P r 2 / \V ffd 2 k a * (y)dy . (3.5) 

The case k = 1 corresponds to Theorem 2.5 in [20] (with A = there). The case k = n 
has been treated in Theorem 3.1 in [29] (see also [28] and Theorem 2.5 in [20]) for any 
a n £ (— § ,0]. We note however that the same proof works for any a n > — So we 
need to consider the intermediate cases k = 2, . . . ,n — 1. 

We deduce (3.5) from the analogous statement for k = n, that is from the following 
inequality 

inf / \f(y)-tf\y\^ dy < C P r 2 f | V/| 2 \y | *£dy , / £ C\B&r)). (3.6) 

£,&IRJB(x,r) " JB{x,r) 

As a consequence of the local representation we have for some a and s = (y, z) 

I d 2 k ^(s)\f{ S )-i\ 2 ds< 
JB(x,r)nn 
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r ra(y)+r—d k (x) 

< C{L) / / \f(y,z)-£\ 2 \a(y)-z\ 2 $dz l dy < 

J \v~ x '\i R r l -k<r Ja) (y)-r-d k (x) 
r rr+d k (x) 

<C / \f(y, a(y) - w) - f| 2 |H^k dwidy , 

J\v-~ x '\m-n,-k< , r Jd k (x)~r 



I2R" 

here we used the following change of variables (y, z) — > (y, w := a(y) — z). Then, since 

i/-H 2 <2(i/-e(y)i 2 + k(y)-eT) , 

where we use the following notation 

. = I^Sf(yMy)-w)\w\^dw 

Id k (x)-r W^dw 
J \y-x \ M n-k<r 

inequality (3.5) follows from estimates (i) and (ii) below. 

(i) We have 

f ( rr+d k (x) o 2a \ 

/ / \f{y, a(y) -w)- £(y)\ Iw]^ dw t dy < 

J \y- x '\ R n-k<r \Jd k {x)-r J 

< C(r + d k {x)) 2 / / \V z f\ 2 \w\nt dw i d y < 

J \y- x '\ M n-k<r \Jd k (x)-r J 

<Cr 2 [ d 2 k a *(s)\Vf\ 2 ds , 

JB(x,r)nn 

here we used the assumption d k {x) < jr as well as inequality (3.6) applied in M k , 
instead of M n , this explains the restriction 2a k > —k. 

(ii) Finally 

rr-\-d k (x) / r \ 

/ / \Z(v)-£\ 2 dy )\w\ 2 ^dw l < 

Jd k {x)-r \J\y-x'\ Rn - k <r j 

<Cr 2 |v,/ + £^-vX(y)| dy\w\% dun < 

Jd k [x)—r J\v—x'\^„_h<r , , UZi 



>d k (x)-r J\y-x>\ Rn _ k <r l=l OZ\ 

<Cr 2 [ dl ak \Vf\ 2 dydz , 

JB(x.r)nn 



/B(x,r)nn 

here we used the standard Poincare inequality on the Euclidean n — k dimensional ball 
of radius r centered at x'. 

The proof of inequality (3.5) is now complete. 

□ 

All the ingredients of the abstract machinery of [24] are now in place. However, 
since bounded domains endowed with the Euclidean metric are not complete manifolds, 
the standard method should be modified as in [20]. In particular we will next prove 
a local weighted Moser inequality as well as a density Theorem which are crucial in 
making the Moser iteration to work in our setting. 

We next prove the following local weighted Moser inequality: 
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Theorem 3.5 (Local weighted Moser inequality) Let n > 2 and ct\ > , at > 

— /or k = 2, ■ ■ ■ ,n. Then, there exist positive constants Cm o-nd r$ such that for any 
u >n + 2A, A := max{ai,ct2, ■ ■ ■ , a n , 0} ; x G ft, < r < r$ and f G Cq°(B(x, r) n ft) 
we have 

r n 

/ \f{y)\ 2 ^ ) \{4^(y)dy< 

JB{x,r)nn " 

<C M r 2 V(x,r)^( [ \Vf\ 2 f[d^(y)dy) (f f f[ d 2 ^{y)dy\ . 

(3.7) 

Proof: Let us first consider the case where d(x) > 7r. Then B(x,r) = B(x,r) C 12 
and for any y G B(x,r) and any k = 1, ■ ■ ■ ,n we have c4(y) ~ dk(x), as in estimate 
(3.2), and the claim follows from the standard Moser inequality, which we recall here: 
There exists a positive constant C such that for any x 6 fi, r > 0, and any v > n if 
n > 3 or any ^ > 2 if n = 2, the following holds true 

2 

[ \f{y)\ 2 ^dy<Cr\-^{f \Vf\ 2 dy)(f fdy\ , 

JB(x,r) \JB(x,r) J \JB(x,r) J 

for all / G Cq°(B(x, r)) (see for example Section 2.1.3 in [30]). Making use of the sharp 
volume estimate in Lemma 3.2, we have 

/ n^(y)l/(y)| 2 ( 1+ ^<Cn^(x)r 2 r-^f/ \Vf\ 2 dy) (f fdy) 

JB{x,r) fj[ £=1 \^B(x,r) / \J B(x,r) J 

<c(f[d^(x)\ 1 Vr-* f / n^(y)iv/i 2 ^ f/ n^(y)/ 2 ^V 

\i=i / \- / ^(*.»-) i= i / V- 75 ^) <=i / 

= C M r 2 V(x,r)^( f f[d 2a >(y)\Vf\ 2 dy)( f f{d 2 ^(y)f 2 dy\ , 

and (3.7) has been proved in case d(x) > 77*. 

Let us now consider the case where d(x) < 7r. Arguing as in Lemma 3.2 this 
corresponds to consider the case where dk{x) < ^r for some k = 1, . . . ,n. In view of 
(3.1) it is enough to prove 

/ \f(y)\ 2 ^d 2 k ^(y)dy< 
JB(x,r)nn 

2 

< CMr'r- 2 -^ 1 ( [ \Vf\ 2 dl a *(y)dy) ( [ f 2 d 2 k ^{y)dy\ " . (3.8) 

In the argument that follows we omit the integral set which is always taken as B(x, r)nf2 
and we define d\x := d 2 k ak (y)dy. First of all making use twice of Holder inequality for 
any v > n + 2a\ , we have 

„+2 "+ 2 "fc 



/ / 2 ( 1+ ^W< J f V ^J dfl J (j d M ) ~ " " +2 < +2 , (3.9) 
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as well as 



(3.10) 

multiplying both sides of inequalities (3.9) and (3.10) we deduce that 



2(n+2a+-i/) 



/ / 2 ^W < (// v ] (/ / 2 ^)^^ r (/ ^ 

Holder inequality also implies that 



(3.11) 



n+2a7-2 



2(n + 2a+) 



J f \ dii<\j f n+2 <- 2 d/j, fdi^j . (3.12) 



To continue we will use the following local weighted Sobolev inequality 

n + 2a+-2 

/ 

J 

V 



2(n + 2ap \ n + 2a + 

I dl a *(y)\f(yT +2 <- 2 dy\ < C s f (y)|V f\ 2 dy . (3.13) 

JB(x,r)nn I JB(x,r)nn 



Then from (3.11), (3.12), (3.13) and (3.3) we get the desired Moser inequality (3.8) 
with 



C M = C s r " +2a t , 

where := max{0, — a&}. It remains to prove (3.13) with a positive constant Cs 
independent of x and r for < r < re,. This is a consequence of (2.28) and the local 
weighted Poincare inequality (3.4) (since for functions / € C^'{B{x,r) n ft) one can 
take £ = in it). It follows that Cm is automatically independent of r if > 0. Hence 
Theorem 3.5 has been proved in this case. 

It remains to show that Cm is independent of r also in the case at < 0. In fact in 
such a case instead of (3.13) we have an even better estimate (see the definition of 
given in (2.9) for q = and n > 3) 

/ dr 2 (y)\f(y)\^dy) <C 8 I ^ d 2 k ^(y)\Vf\ 2 dy, 

for some positive constant tq and for any x <E Q, < r < vq, f <E C^°(B(x, r) fl fl), with 
a positive constant C5 independent from x and r. Whence 

n-2 4a k 

([, , 4 Qfc (y)l/(y)l^^) " <Cs( max 4(y)) " / ^ d 2 k ^(y)\Vf\ 2 dy< 

\JB{x,r)nn I \y£B{x,r)nn I JB(x,r)nn 
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<C s r— d 2 h a *(y)\Vf\ 2 dy , 

JB(x,r)nn 

since d k (y) < dk(x) + r < (1 + r y)r. If we use the above inequality in the place of (3.13) 
we get Cm = Cs, that is Cm is independent of r. 
The proof of Theorem 3.5 is now complete. 

□ 

Finally let us prove the following density Theorem, which as explained in [20] is 
crucial for the Moser iteration to work on bounded domains. Let 

H 1 ^; cf>l) = {ve L 2 {Q,(t)ldx) : f (f>l\Vv\ 2 dx < +oo}, 

Jn 

with norm defined by |H|^i := Jn<i>i(v 2 + \ Vv\ 2 )dx. 

4>i 

Theorem 3.6 (Density theorem) Let n > 2. Suppose that <pi satisfies 

Cl d^(x)...d^(x)<Mx)<c 2 d^(x)...d^(x), 
for with c\, c 2 positive constants and a k > — k = 1, . . . ,n. Then, the 

C£°(Q) functions are dense in iJ 1 (f2; cf>\). 

Proof: The special case ct k = for k = 2, . . . , n was treated in Theorem 2.11 of 
[20]. First of all from Theorem 7.2 in [25] it is known that the set C°°(f2) is dense in 
H 1 (VL; 4>\) . Thus for any v G H 1 ^;^ 2 ) there exists v m £ C°°(fi) such that for any 
e > we have \ \v — v m \\ H i < e if m > m(e). Let us choose w := v m i e \ and let us define, 
for any k = 1, ■ ■ ■ , n and any j > 1, the following function 



V k {x) 



if d k (x) < , 

^ 1 if d k {x) > \ . 

Then wi :=w]Jt=i 4>{ G C^(fl), and 

n „ n 

\\w-uS\\ Hl =|Kl-n^)llm <2 (w 2 + \Vw\ 2 )(l-\l^ k ) 2 ( t ) 2 (y) dy+ 
1 k=i 1 Jn k=i 

™ 2 £lwiiV?(y) d v< 
1 k=i 

<2/ {w 2 + \Vw\ 2 )cp 2 dy + 2Y / ,9 » ,\^ dl ak (y)dy . 



Now as j — > cx) it is clear that the first term in the right hand side goes to zero 
since w 6 H^. We next show that also the second term goes to zero. Recalling that 
|V4| < 1 

/ J^%d 2 k ak (y)dy < c l ^f^ I t 2 ^- 2 t^dt < 

Jji<d k (v)<i d 2 (y)(ln(j)y k (ln(j)) 2 J^<t<l 
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- (l n (j 2 j2a k -2+k Jj^ <t< , ~ l n(j ) 

as j — > oo for any 2a/c — 2 + k > 0, and this completes the proof. 

□ 

At this point we have all the ingredients needed in order to apply Moser iteration 
technique up to the boundary, as adapted on bounded domains in [20], to the operator 



1 n 8 ( n f)v \ 

C ^--=- r,n E 7^ iK^fc )+Xiv , (3.14) 



or equivalently to the degenerate elliptic operator C^, defined in (2.37). 
In fact one can prove the following result 

Theorem 3.7 (Parabolic Harnack inequality up to the boundary) Let n > 2, 

ft C M n , be a smooth bounded domain, Ai S M and > —^^ L , for k = l,---,n. 
Then there exist positive constants Ch and R = R(Q) such that for i 6 fi, < r < R 
and for any positive solution v(y, t) of 

— = -C a v in {B(x, r)nn}x (0, r 2 ) , (3.15) 
the following estimate holds true 

ess ™P Me{B{x r _ )nn}x{ ^ ^ )V (y,t) < C H ess 'mf Me{B{Xt r )m}x( s^ ra) v(y,t) . 

Here we use the following definition of solutions: 
Definition 3.8 By a solution v(y,t) to (3.15), we mean a function 

v € C 1 ({0, r 2 ); I? (b(x, r)nQ,f[ df k (y)dy^j j nC° ^(0, r 2 ); H 1 (b(x, r) fl fi, f[ af k (y)dy 
such that for any $ G C°((0, r 2 ); Cq°(B(x, r) n fi)) and any < ii < t 2 < r 2 we have 



pf 2 r \ n r) /)(T) I n 



Let us note that Theorem 3.7 is sharp, in the sense that the same statement does 
not hold true if at < — for some k = 1, . . . , n as explained in [20]. 

The parabolic Harnack inequality up to the boundary for the Schrddinger type 
operator C defined in (1.9) and stated in Theorem 1.2 is proved as follows: 

Proof of Theorem 1.2: Clearly Theorem 3.7 applies also to the operator instead of 
C a . Hence Theorem 1.2 is a consequence of Theorem 3.7 for and of the change of 
variables v = ucj)^ 1 see (2.36). 

□ 
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^,From the parabolic Harnack inequality in Theorem 3.7 we deduce, as in [20], the 
sharp two-sided estimates for the heat kernel l a associated to the elliptic operator C a 
defined in (3.14) under Dirichlet boundary conditions. That is 

„ n 

v(x,t) := / l a (t,x,y)v (y)Y[d 2 i a '(y)dy 
Jn t=i 

satisfies vt = —C a v in (0, oo) x Q, v(0, x) = vq(x) on Q and v = on (0, oo) x d\£l. We 
then have 

Theorem 3.9 Let n > 2, £1 C lR n , be a smooth bounded domain, Ai G 1R and at > 
f or k = 1, ■ ■ ■ ,n. Then there exist positive constants C\, C2, with C\ < C2, and 
T > depending on 17 such that 

n 1 - 1 2 

CiHidiix) + Vt)- ai (di(y) + Vty^t-^e- 02 ^- <l a (t,x,y) < 

i=l 

n I — I 2 

< c 2 H(di(x) + Vi)-*^) + Vt)- a n-^e- c ^^- 

1=1 

for all x,y G £1 and < t < T. 

Finally from the global upper bound in (2.38), arguing as in Theorem 6 of [7] (see 
also Proposition 4 in [8] as well as the proof of Theorem 1.2 in [20]), one can deduce 
an analogous lower bound for large times, thus obtaining the following sharp long-time 
asymptotics of the heat kernel 

Theorem 3.10 Let n > 2, U C M n , be a smooth bounded domain, Ai G M and 
ctk > — for k = 1, • ■ ■ , n. Then there exist positive constants C±, C2, with C\ < C2, 
and T > depending on 17 such that 

Cie- Alt < l a (t,x,y) < C 2 e- Xlt 

for all x, y G Q, and t >T. 

^From Theorems 3.9 and 3.10, making use of the equivalence (2.39) as well as of 
assumption (1.6), we get the corresponding result for the Schrodinger operator C stated 
in Theorems 1.3 in the Introduction. We omit further details. 

As we have already mentioned integrating the sharp two-sided estimates for h(t, x, y) 
in Theorem 1.3 with respect to the time variable, one can deduced estimates on the 
Green function for the Schrodinger operator C defined in (1.9) in the case Ai > 0. Some 
explicit examples of sharp two sided Green function estimates are given in Theorem 
4.11 in [20]. 

4 Applications 

In this Section we give some examples of singular potentials V for which the results 
of the present work apply; that is, we give examples of potentials V for which the 
generalized first eigenvalue is not —00, and the corresponding first eigenfunction is 
bounded from above and below uniformly by some power of the distance function. We 
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should stress that the asymptotics of cf>i for the examples that follow is a consequence 
only of the maximum principle as used in [5]. We will present the detailed argument 
for example III in the Appendix; the other cases can be treated similarly. 

To this end let us first prove that the sum of two potentials having disjoint sin- 
gularity sets and finite generalized first eigenvalues, also has finite generalized first 
eigenvalue. 



Lemma 4.1 Let V; L , i = 1,2 be such that Vi 6 Lj oc (Q) n Lf£ c (£l \ Si), where Si are 
compact subsets of 17 such that Si n S2 = and 

«ecg°(n) J n u 2 dx 
Assuming that Ai(V^) > —00 for i = 1,2 then Xi(Vi + V2) > — 00. 

Proof: Let us take (p 6 C°°(fi), < 92 < 1, such that </? = 1 in finfii and (p = in f2\f2i 
where S7i is a neighborhood of S±, that is f^i := {2; G f2 : dist(x, Si) < (5} for some small 
5 > 0, and Oi is a slightly bigger neighborhood of Si, thus S\ C f2i C Q\. Whence 
SI2 := is a neighborhood of S2 and ^2 := is a slightly bigger neighborhood 

of S*2, thus O2 C ^2- Whence for any u £ Cq°(Q) we have u = mp + u(l — ip) =: u\+U2- 
By elementary calculations we have that 



\Vu\ 2 dx = / \V(uip + (1 - (p)u)\ 2 dx > / \Vui\ 2 dx+ / \Vu 2 \ 2 dx - K u 2 dx 
Jn Jn Jn Jn 

for a suitable positive constant K. Then for V := V± + V2 we have 

j (|Vu| 2 - VV) dx > 

> (|Vm| 2 - V lU \ + |Vu 2 | 2 - V 2 u 2 2 + u 2 {V W 2 + V 2 (l -^) 2 -V- if)) dx = 
= (\V Ul \ 2 - V lU 2 ) dx + (\Vu 2 \ 2 - V 2 v%) dx+ 



dx > 



+ f u 2 W^ 2 - 1) + V 2 ((l - <f) 2 - 1) - K 
Jn 1 J 

> (Ai(Vi) + Ai(^ 2 )- ll^ill^^-H^I^)-^) / u 2 dx 



□ 



We now present some concrete examples. 

Example I Our first example is motivated by [15], [16], [17] and deals with multipolar 
inverse-square potentials. Let n > 3, f2 C M n be a smooth bounded domain from which 
we have removed m points xi,..., x m and 

m 

V(x) = £ 



2 

i=i 



for < Ci < — P . We note that differently from [15], we may take in each one of the 
inverse-square potentials the critical Hardy constant. This is due to the fact that we 
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study the Schrodinger operator —A — V on a bounded domain. In such a case one can 
prove that 0i(x) ~ Il^=i \ x ~ Xi\ l3t dist(x,diQ) with 



_ 2 - n + y/(n - 2f - 4g 
Pi-- 2 

see Lemma 7 in [13] and Theorem 7.1 in [11] on one hand and the elliptic regularity on 
the other. In fact the function f(x) := \x — Xi\^ satisfies the equation A/+ ^rB^z f = 
in M n \ {xi}. We only need to check that Ai > — oo. This follows from Lemma 4.1, 
which clearly can be generalized to a finite sum of potentials, and the improved L 2 
inequality given in [31] for a single inverse-square potential. Consequently Theorems 
1.2, 1.3 and 2.1 apply with r*. = for k = 2, . . . , n — 1 and T n = {x\, . . . ,x m }. Note 
that in this case d" n (x) stands for n^Li \x— Xj|^ 1 , that is, we have m different sets of the 
same codimension n in the boundary of f2, where <f>i may present different degeneracies. 

Example II Let n > 4, Q = Br \ E, for some R > 1, where 

E := {x e TR n : x\ + x\ = 1, x 3 = . . . = x n = 0} 
B R := {x e M n : x\ + . . . + x 2 n < R 2 } and 



4disr(x,5S R ) 4 disr(x,£?) 



3 — n 



In such a case one can easily prove that 4>\{x) ~ dista (x, 8Br) dist^ - (x, E 1 ), see [13]. 
The fact that Ai > — oo follows making use of Lemma 4.1, from the improved L 2 
inequality given in [4] for the inverse-square distance to OBr and the one given in [13] 
or [19] for the inverse-square distance to the set E having codimension n — 1. Theorems 
1.2, 1.3 and 2.1 now apply to —A — V with a\ = 1/2 and a Tt _i = (3 — n)/2, whereas 
all the other a^s are zero. 

Example III Let n > 2, and Q C M n be a bounded domain such that dQ = <9i$7, 
that is, the boundary of Q has codimension one. We now take 

V(x) ; ; 



4dist^(x,^) 

By the results of [4] we have that Ai > — oo under appropriate regularity assumptions 
on dd. We recall also that 4>i(x) ~ dista (x, d£l), as shown in [13]; see also Appendix, 
where we will provide a self-contained proof based only on the maximum principle. 
Therefore Theorems 1.2, 1.3 and 2.1 apply to —A — V with oi\ = 1/2, whereas all the 
other afc's are zero. We note that this improves the corresponding Theorems in [FMT1] 
removing the convexity assumption, under which it is known that Ai > (see [4]). 

Example IV Let n > 3, ffl C M n be a smooth bounded domain containing the origin, 
f} = jy\{0}, and 

1 1 (n-2) 2 1 

{x) ~ 4dist 2 ( x ,d 1 n) + 4 |x| 2 • 

The fact that Ai > — oo may be deduced from Lemma 4.1 making use of the L 2 improved 
Hardy inequality in [31] for the inverse-square potential l/|x| 2 and of the one in [4] for 
the inverse-square potential involving the distance to di$l. In this example we have 

1 2 — n 

(j>i(x) ~ distz (x, diO,)\x\^~ . Whence Theorems 1.2, 1.3 and 2.1 apply to — A — V with 
cci = 1/2 and a n = (2 — n)/2, whereas all the other a^'s are zero. 
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Example V Let n > 3, and B\ C M n be the unit ball. For — < a < we consider 
the operator C = — L — V where 

8 ( 8U \ 1 .O-n,* p X 

ij"=l 1 \ V 

and V(x) = — a ^p~ 2 ' > • The operator L is easily seen to be uniformly elliptic and in 
fact 

On the other hand if 

du du 



I f J2 M a 



using the change of variables u = \x\ a v a straightforward calculation shows that 

J Bi ^ E ai ^§^§^ ~ Vu ^j dx = Jg^W 2 " + \\ x \ a+2{ ± ~ 8ij))v Xi v Xi dx. 

Using (4.1) one can easily see that Ai > 0. 

In this example we have that <p\(x) ~ dist(x, dB\)\x\ a . We can apply our results 
with a± = 1 and a n = a, whereas all the other q^'s are zero. 

5 Appendix 



In this Appendix we consider the operator C := —A — ^-jt^ which corresponds to 
Example III of Section 4 and we will prove that the corresponding eigenfunction <\>\ is 
such that 4>\(x) ~ dz(x), where d(x) := dist(x, <9f2). 

S'tep /: Existence of <p\ in a suitable energy space. Let 77 £ C 2 (0) be a function such 
that 77(3;) = d 1 / 2 (x) near the boundary, say, d(x) < eo, and r/(x) > Co > for > eo- 
Let Hj(£i) be the closure of CQ°(f2) functions under the norm 



|^ ■- J^d(\Vv\ 2 + v 2 )dx . 



This norm is equivalent to the norm \\v\\ := J n r] 2 (|Vw| 2 + v 2 ) dx. Changing variables 
by u = rjv, in 



U\Vu\ 2 -f^)dx 



-00 < Ai = inf 



«ecg°(n) jf^ 2 ^ 
we get the equivalent inequality 

. f ^(^|V^| 2 -(^r/ + gy)dx 
-00 < Ai = inf — — . (5.1) 

vec ~(n) fnV 2 v 2 dx 

2 

Using the fact that r] Ar] + 33 G L°°(fi) as well as the following estimate 

[ ^l v ^dx < C [ d{\Vv\ 2 + v 2 )dx , 1; G c °°(n), 
Jn d Jn 
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which was established in Proposition 5.1 in [21], we find that for every e > there 
exists an M e > such that for all v G C§°(Q), 

[ (r]Ar] + ^-)v 2 dx <e [ ij 2 \Vv\ 2 dx + M e [ rj 2 v 2 dx. (5.2) 
Jn 4a Jn Jn 

In the sequel we will establish the existence of a function tpi G H^(Q) which realizes 
the infimum in (5.1). To this end let w/. be a minimizing sequence normalized by 
In V 2w k^ x = 1- Then, using (5.2) we can easily obtain that the sequence Wk is bounded 
in H\. Therefore there exists a subsequence still denoted by wu such that it converges 
H\ - weakly to ipi, and in addition we have the following strong convergence, for 
k — ► OO, 



and 



(rjAri + —)wldx ^ j^(r)Ar) + —)ipjdx, 

j 2 2 / 2 2 

/ r) w k dx — > / 77 ipidx. 
Jn Jn 



Using the lower semicontinuity of the gradient term in the numerator of (5.1) the result 
follows. 

Step II: An auxiliary estimate. For 8 > small enough we set 

L (d\Vv\ 2 - Ad 4)dx 

Mi(^):= inf Us V - - _ , 2 > , (5.3) 
vec~(n s ) In s dv 2 dx 

where tt$ = {x G SI s.t. dist (x, 9f2) < <5}. 
We will show that 

Hi(Cls) —> +00, as <5 — ► 0. (5.4) 

Our starting point is the inequality 

j Qs (|Vu| 2 - l^j dx>\ ^-u 2 dx , u € C^(Q 5 ), (5.5) 

for any < 5 < So, for some 8q small enough, where X{t) := (1 — hit) -1 . To prove this 
one starts with the obvious relation 



~ Jn, \2d 2d J 



dx . 



Expanding the square, integrating by parts and using the fact that \dAd\ can be made 
arbitrarily small in Qg, for 8 sufficiently small, the result follows. Changing variable as 
usual by u = d?v, inequality (5.5) is equivalent to 

/ [ d\Vv\ 2 - Ad—) dx>-[ ^-^-v 2 dx , v e C^(Q 5 ) • (5.6) 
Jn 6 \ 2 J 8 Jn s d 



Therefore 



\X\8) X 2 (5)In s (d\Vv\ 2 -Ad v 4)dx (d\Vv\ 2 - Ad^) dx 

8 8* ~ 8* I^^vHx ~ In^ddx 
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from which (5.4) follows. 

Step III: Asymptotics of (f>\. The lower bound C\d l l 2 (x) < <j)\(x) is a consequence of 
the maximum principle and is derived in Lemma 7 in [13]. 

We will obtain the upper bound using maximum principle in a suitably small neigh- 
borhood of the boundary. Let^i(x) = <pi(x)/d 1 ^ 2 . Then, for Ev := —div(dVv) — ^v — 
Xidv, we have that Eipi = 0. Moreover, we have that 

E(l - Cd) = C - ^ + d + ^-Ad + CXn^j > 0, 

in Qg for 5 small enough and C > big enough. We next choose f3 > big enough so 
that 

Mx) < - Cd) on dfls- 

Let g(x) := 4>i{x) — — Cd) and g + := max{0,g}. We clearly have that 

/ 9 + Eg<0, 
Jn s 

from which it follows that 

fn s (d\Vg + \ 2 - ^(g+) 2 )dx 

J ns d(g+) 2 dx ~ 

This contradicts (5.4) unless g + = from which it follows that <f>(x) < (3d l l 2 (x). 
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